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1.  INTRODUCTION 

Throughout  X  =  {X(t),  t  6  T}  is  a  real  symmetric  o-stable  (SaS)  process  with 
0  <  a  <  2  and  T  an  interval  on  the  real  line;  i.e.  all  finite  linear  combinations 

,  a  X(t  )  are  SaS  random  variables. 
n  =  i  n  tr 

X  is  Markov  if  for  all  s  <  t  in  T,  the  conditional  distribution  of  X(t)  given  X(u), 
u  <  s,  "coincides"  with  the  conditional  distribution  of  X(t)  given  X(s)  alone,  in  the  sense 
that  for  ?ny  u-  <  s,  i=l,...,n,  and  any  Borel  set  E  the  equality  P{X(t)  €  E|X(u1), 

X(u  ),  X(s)}  =  P{X(t)  e  E  |  X(s) }  holds  with  probability  1.  As  is  well-known  the 
Markovian  property  is  equivalent  to  the  conditional  independence  of  the  past  <r{X(u), 
u  <  t}  and  future  cr{X(u),  u  >  t}  cr-fields  given  the  present  <r{X(t)},  and  thus  it  is 
symmetric  in  time  and  could  be  defined  by  requiring  that  for  all  t  <  s  in  T  the  conditional 
distribution  of  X(t)  given  X(u),  u  >  s,  "coincides"  with  the  conditional  distribution  of  X(t) 
given  X(s)  alone.  Conditional  distributions  of  non-Gaussian  stable  processes  are  generally 
very  difficult  to  compute  (and  generally  not  stable)  and  it  is  thus  not  easy  to  check  for  the 
Markovian  property.  For  this  reason  we  introduce  a  weaker  Markovian  property  which  is 
amenable  to  some  analysis  and  which  concentrates  on  regressions. 

For  1  <  a  <  2  we  have  <S\ X(t)  |  <  oo  and  we  say  that  X  is  left  weakly  Markov  if  for  all 
s  <  t  in  T  with  probability  1, 

<?{X(t)|X(u),u<s}  =  *{X(t)|X(s)}, 

and  right  weakly  Markov  if  for  all  s  <  t  with  probability  1, 

£{X(s)  | X(u),  u  >  t}  =  £{X(s) | X(t)}. 

In  the  Gaussian  case  o=2  the  left  and  right  weak  Markovian  properties  are 
equivalent,  and  they  are  also  equivalent  to  the  Markovian  property.  Furthermore,  there  is 
only  one  stationary  Gaussian  process  which  is  Markov,  namely  the  Ornstein-Uhlenbeck 
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process  with  covariance  function  R(t)  =  R(0)e  In  contrast  there  are  stationary 

non-Gaussian  S aS  processes  with  1  <  a  <  2  which: 

(i)  are  left  weakly  Markov,  without  being  right  weakly  Markov,  and  vice  versa 
(cf.  Section  7); 

(ii)  are  left  and  right  weakly  Markov  without  being  Markov,  e.g.  the  sub- 
Ornstein-Uhlenbeck  processes  {Corollary  4.2); 

(iii)  are  Markov,  namely  the  SaS  Ornstein-Uhlenbeck  processes  in  (2.14)  whose 
covariation  function  is  the  nonsymmetric  double  exponential  function  in 
(2.15); 

(iv)  have  the  symmetric  covariation  function  R(t)  =  R(0)e-^^^  but  are  neither 
left  nor  right  weakly  Markov,  namely  the  harmonizable  process  in  (6.3). 

T»vo  distinct  S aS  stationary  Markov  processes  are  identified  in  this  paper.  These  are 
the  right  and  the  left  S aS  Ornstein-Uhlenbeck  processes,  which  can  be  represented 
respectively  as  decreasing  and  increasing  time  changes  of  SaS  Levy  motion  (cf.  (2.14)  and 
(3.3)),  or  as  nonanticipating  and  fully  anticipating  moving  averages  of  SaS  Levy  motion 
(Theorem  5.1  and  (5.4)),  and  are  the  stationary  solutions  of  certain  first  order  stochastic 
differential  equations  driven  by  SaS  white  noise.  Even  though  there  might  be  further  SaS 
stationary  Markov  processes,  none  is  currently  known.  Such  processes  are  not  sub- 
Gaussian  (Corollary  4.2)  or  harmonizable  (Theorem  6.1);  and  they  are  neither 
nonanticipating  nor  fully  anticipatory  invertible  moving  averages,  as  the  left  and  the  right 
SaS  Ornstein-Uhlenbeck  processes  are  the  only  such  SaS  moving  averages  (Theorem  5.1 
and  page  33).  Finally,  neither  one  of  their  pairwise  conditional  distributions  can  be 
o-stable  and  symmetric,  since  the  left  and  the  right  SaS  Ornstein-Uhlenbeck  processes 
are  again  the  only  ones  possessing  this  property  (Theorem  3.1).  This  may  explain  the 
difficulties  in  constructing  other  SaS  stationary  Markov  processes,  if  indeed  there  are  any. 

Without  requiring  stationarity,  the  Gaussian  case  is  still  quite  simple:  All  Gaussian 
Markov  processes  are  essentially  time  changes  of  Brownian  motion,  see  Tismoszyk  (1974), 
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Borisov  (1982)  and  Wong  and  Hajek  (1985).  For  non-Gaussian  SaS  processes  with 
1  <  a  <  2  the  picture  is  more  complex  and  rich.  A  necessary  condition  for  left  weak 
Markovianness  is  given  in  Theorem  2.1  in  terms  of  the  covariation  function,  and  its 
solution  is  found,  i.e.  in  the  nonstationary  case  (2.5)  and  its  generalization,  and  in  the 
stationary  case  (2.9).  While  all  time  changes  of  Levy  motion  have  covariation  function  of 
this  form  (i.e.  (2.5)),  and  are  in  fact  Markov,  they  do  not  exhaust  the  class  of  SaS  processes 
with  covariation  function  of  this  form,  e.g.  Levy  bridge  (see  Example  2.1). 

Time  changes  of  Levy  motion  are  considered  in  Section  3  where  they  are  shown  to  be 
the  only  SaS  Markov  processes  whose  pairwise  conditional  distributions  are  stable  and 
symmetric  (Theorem  3.1).  In  the  non-Gaussian  stable  case  there  is  also  a  marked 
asymmetry:  The  SaS  Markov  processes  whose  right  to  left  and  left  to  right  pairwise 
conditional  distributions  are  stable  and  symmetric  are  few  and  trivial  when  1  <  a  <  2, 
whereas  in  the  Gaussian  case  a=  2  they  coincide  with  the  entire  class  of  Gaussian  Markov 
processes. 

An  auxiliary  result  of  independent  interest  is  given  in  Proposition  3.1  and  Corollary 
3.1,  characterizing  the  stability  of  the  conditional  distribution(s)  of  random  variables  that 
are  jointly  SaS. 

Sub-Gaussian  processes  are  left  (right)  weakly  Markov  if  and  only  if  they  are 
essentially  time  changes  of  sub-Brownian  motion,  except  for  trivial  cases,  and  they  are  not 
Markov  (Theorem  4.2).  In  particular,  the  only  weakly  Markov  stationary  SaS 
sub-Gaussian  processes  are  the  sub-Ornstein-  Uhlenbeck  processes  (Corollary  4.2). 

Sections  5  and  6  consider  two  specific  classes  of  stationary  SaS  processes,  moving 
averages  and  harmonizable  SaS  processes.  It  is  shown  that  in  the  case  of  either 
nonanticipating  or  fully  anticipatory  invertible  moving  averages,  the  weak  Markov 
property  cannot  exist  without  full  Markovianness  and  it  is  realized  only  by  the  right  and 
the  left  SaS  Ornstein-Uhlen beck  processes  correspondingly  (Theorem  5.1  and  page  33). 
In  sharp  contrast  to  the  Gaussian  case  o=2,  it  turns  out  that  for  the  stable  case 
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2.  GENERALITIES 


A  SaS  process  X  can  be  represented  by  an  integral  of  the  form 


X(t)  =  f ( t , u )  dZ(u),  t  G  T, 

Jit 


where  Z  is  a  SaS  random  measure  on  some  <7-finite  measure  space  (U,E,//)  fi.e.  Z  is  an 
independently  scattered  cr-additive  set  function  on  E  —  {Eg  E,  //(E)  <  oc}  and 

ji 

£exp{irZ(E)}  =  exp{-//(E)|r|  a}  for  E  G  E^)  and  {f(t, - ),  t  G  T}  C  Lft(U,E,//)  (Kanter 
(1972),  Kuelbs  (1973)  and  Hardin  (1982)).  All  quantities  are  real-valued  (except  in 
Section  6,  where  complex-valued  processes  are  discussed)  and  //  is  called  the  control 
measure  of  Z.  For  every  g  G  Lq(//)  the  integral  / gdZ  is  a  SaS  r.v.  with  £exp{ir/gdZ}  = 
exp{-|r|°/|g|ad/'}  and  is  linear  in  g.  Specific  examples  of  such  representations  of  SaS 
processes  will  be  considered  in  Sections  5  and  6. 

The  covariation  function  of  X  is 

(2.2)  R(t,s)  =  Cov[X(t),  X(s)j  =  y  f(t,u)  f(s,u)<a_1>d//(u), 

where  x<ci>=  ( x  |  ^  sgn(x),  and  it  does  not  depend  on  the  specific  representation  of  X  (for 
more  on  the  covariation  see  Section  3).  In  the  Gaussian  case,  a=2,  the  covariation  is  a 
multiple  of  the  covariance,  R(t,t)  determines  the  distribution  of  X(t),  and  the  numbers 
R(t,t),  R(s,s),  R(t,s)  determine  the  joint  distribution  of  X(t),  X(s);  thus  knowledge  of  R  on 
T*T  determines  the  distribution  of  the  (zero-mean)  Gaussian  process  X.  In  the 
non-Gaussian  SaS  case  1  <  a  <  2,  the  covariation  is  not  generally  a  symmetric  function  of 
its  arguments  and  is  linear  only  in  the  first  argument,  R(t,t)  determines  the  distribution  of 
X(t),  but  the  numbers  R(t,t),  R(s,s),  R(t,s),  R(s,t)  do  not  generally  determine  the  joint 
distribution  of  X(t),  X(s).  Thus  knowledge  of  the  covariation  function  R  on  TxT  generally 
does  not  determine  the  bivariate  distributions  of  the  SaS  process  X.  Still,  as  we  shall  see, 
the  covariation  function  plays  a  role  partially  analogous  to  the  role  played  by  the 
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covariance  function  in  the  Gaussian  case. 

A  basic  result  on  the  left  weak  Markovian  property  is  the  following: 

Theorem  2. 1 .  X  is  left  weakly  Markov  if  and  only  if 

(2  3)  Cov[X(t)-{$if}x(s),  Y|  =  0 

for  all  s  <  t  and  all  Y  6  sp{X(u),  u  <  s},  where  the  closure  is  in  probability.  If  X  is  left 
weakly  Markov  then 

(2.4)  R(t3,t2)  R(Mi)  =  R-(t3’ti)  R(t2’t2)  for  all  t1  <  t2  <  tg. 

A  covariation  function  R  with  R(t,s)  f  0  for  alls  <  t  in  T  satisfies  (2.4)  if  and  only  if  it  is  of 
the  form 

(2.5)  R(t,s)  =  H(t)  K(s)<q“1>  for  all  s  <  t, 

where  the  functions  K,  H  are  unique  up  to  a  multiplicative  constant ,  have  the  same  sign  and 
K(t)/H(t)  is  positive  and  nondecreasing  on  T. 

Proof.  It  is  known  (see  Kanter  (1972))  that 

*{x(t)|x(s))  =  Rjf|)x(S). 

Therefore  X  is  left  weakly  Markov  iff 

<T{X(t)|X(u),u<s}  =^||x(s),  V  s  <  t, 

and  by  [3,  Proposition  1.5],  a  necessary  and  sufficient  condition  for  this  is  (2.3)  for  all  s  <  t 
and  Y  e  sp{X(u),  u  <  s}. 

Now  if  X  is  left  weakly  Markov,  taking  Y  =  X(u),  u  <  s,  we  obtain 


R(t,u)  =  R(s,u),  V  u  <  s  <  t, 


which  is  (2.4).  The  general  form  (2.5)  of  the  solution  of  (2.4)  is  obtained  as  in  Borisov 
(1982)  by  taking,  for  some  interior  point  tQ  of  T,  K(t)<Qi-1>=  R(tQ,t)  for  t  <  t  , 

=  R(t,t)  R(tQ,t0)/R(t,t0)  for  t  >  tQ,  and  H(t)  =  R(t,t)/R(tQ,t)  for  t  <  tQ, 

=  R(t,tQ)/R(tn,t0)  for  t  >  tQ.  Since  by  (2.2),  R(t,t)  >  0  and  by  assumption  R(t,t)  f  0,  it 
follows  from  0  <  R(t,t)  =  H(t)K(t)<Q^>  that  K  and  H  have  the  same  sign  at  each  point. 
Also  from  (2.2)  and  Holder's  inequality  we  obtain 

(2.6)  | R(t,s) |  <{R(t,t)}1/a{R(s,s)}1_1/a 

and  substituting  from  (2.5)  we  have  |  K(s)/H(s)  |  <  |  K(t)/H(t)  | .  Since  KH-1  is  positive,  it 
is  nondecreasing  on  T.  Conversely,  (2.5)  implies  (2.4)  immediately  and  the  property  KH_I: 
nondecreasing,  is  needed  to  show  that  R  given  by  (2.5)  is  covariation  function.  The 
simplest  way  of  showing  this  is  by  constructing  a  SaS  process  with  covariation  (2.5),  as  was 
done  in  the  Gaussian  case  in  Wong  and  Hajek  (1985),  p.  64.  Indeed,  using  the  time  change 
r(t)  =  {K(t)H  ^(t)}0^  (nondecreasing),  and  the  SaS  Levy  motion  L  =  (L(t),  t  >  0} 
which  has  stationary  independent  increments,  L(0)  =  0,  and  £exp{ir[L(t)-L(s)]}  = 
exp{-|  r|  a|  t-s  | },  we  can  introduce  the  SoS  process 

(2.7)  X(t)  =  H(t)L(r(t)) 
whose  covariation  function  is  for  s  <  t, 

Cov[X(t),X(s)J  =  H(t)H(s)<^>  Cov[L(r(t)),  L(r(s))J 

=  H(t)H(s)<c^1>  r(s)  =  H(t)H(s)<Q>_1> 

(2.8)  =  H(t)K(s)<Q^1>=  R(t,s) 

since  K(t)H(t.)  >0.  □ 


In  the  Gaussian  case  a=2,  the  covariation  is  linear  in  its  second  argument  (as  well  as 


in  its  first),  and  the  necessary  condition  (2.4)  is  also  sufficient;  thus  when  R(t,t)  t  0,  t  e  T, 
conditions  (2.3),  (2.4),  (2.5)  and  (2.7)  are  ail  equivalent,  and  all  Gaussian  Markov 
processes  are  time  changes  of  Brownian  motion.  However,  in  the  non-Gaussian  SaS  case 
with  1  <  a  <  2,  generally  the  covariation  is  not  linear  in  its  second  argument  and  the 
necessary  condition  (2.4)  is  not  sufficient.  Also,  while  the  time  changes  of  SaS  Levy 
motion  (2.7)  have  covariation  function  of  the  form  (2.5),  they  do  not  exhaust  the  class  of 
SaS  processes  with  covariation  function  of  the  form  (2.5). 

/ 

Example  2. 1  Levy  bridge. 

Again  let  L  be  the  Levy  motion,  and  let  B(t)  =  L(t)  -  tL(l),  0  <  t  <  1.  This  is  one  of 
the  possible  generalizations  of  the  Brownian  bridge  to  the  SaS  case,  a<2.  It  is 
starightforward  to  check  that  for  this  process 

( (l-tlsKl-s)0-^  s®-1]  if  0  <  s  <  t  <  1, 

R(t.s)  =  Cov[B(t),  B(s)]  =  !  . 

( t  ( 1— s)  [( 1 — s)  +  s  ]  if  0  <  t  <  s  <  1. 

Moreover,  B(t)  is  easily  seen  to  satisfy  the  condition  (2.3)  for  any  Y  =  ^j=iajX  (m), 
Uj<s,  for  i=l,2,...,k,  and,  therefore,  for  any  Y  G  sp{X(u),  u  <  s}.  This  process  is, 
therefore,  left  weakly  Markov  and,  in  fact,  two-sided  weakly  Markov,  since  its  right  weak 
Markovianness  can  be  established  similarly. 

The  Levy  bridge  B(t)  is  an  example  of  a  two-sided  weakly  Markov  SaS  process 
which  is  not  a  time  changed  Levy  motion  (see  (2.12)).  Other  examples  of  such  processes 
are  furnished  by  the  sub-Gaussian  SaS  processes  (see  Section  4). 

The  process  B(t)  is  probably  not  Markov.  It  is  interesting  to  note  that  another 
possible  generalization  of  the  Brownian  bridge,  namely  B'(t)  =  ( 1— t) L[t/(  1— t)] ,  0  <  t  <  1,  is 
clearly  distinct  from  B(t)  when  a  <  2  !  B'  is  a  Markov  process,  and  its  covariation 
function  is  given  by 


WWW 


wv 


R'(t,s)  =  Cov[B'(t),  B'(s)] 


(l-t)s(l-s)  if  0  <  s  <  t  <  1, 

t(l-s)^1  if  0  <  t  <  s  <  1. 


Not  much  seems  to  be  known  about  the  role  of  these  processes  (if  any)  in  the  weak 
convergence  of  empirical  processes.  □ 


The  solution  of  (2.4)  in  the  general  case,  i.e.  without  the  condition  R(t,t)  t  0  on  T, 
can  be  obtained  as  in  the  Gaussian  case  (Timoszyk  (1974),  Borisov  (1982)),  in  the  form 
(2.5)  on  a  finite  or  denumerable  union  of  disjoint  squares  around  the  diagonal  of  T*T  (and 
zero  elsewhere). 

When  X  is  stationary  (T  =  3^)  then  R(t,s)  depends  only  on  t-s  and  we  write  R(t,s) 
=  R(t-s).  When  o=  2  the  converse  is  also  true,  but  this  is  not  generally  true  when 
1  <  a  <  2.  When  R(t,s)  =  R(t-s)  for  all  t,s  e  IR1,  we  say  that  X  is  covariation  stationary. 

In  the  presence  of  stationarity  Theorem  2.1  reduces  to  the  following  simpler  form. 

Corollary  2.1 .  Let  T  =  IR*.  If  X  is  covariation  stationary  and  left  weakly  Markov ,  then 
for  some  0  <  A  <  oo, 

(2.9)  R(t)  =  R(0)e~At  for  all  t  >  0. 


IfX  is  stationary ,  then  it  is  left  weakly  Markov  if  and  only  if 
(2.10)  Cov(X(t)  -  e~AtX(0),  Y]  =  0 

for  some  0  <  A  <  oo  and  all  t  >  0,  Y  6  sp{X(u),  u  <  0}. 

Proof.  If  X  is  covariation  stationary  and  left  weakly  Markov,  then  (2.4)  is  satisfied  and 
can  be  written  in  the  form 


R(u)R(v)  =  R(u+v)R(0)  for  all  u,v  >  0. 


Since  by  (2.6),  |  R(t)  |  <  R(0),  V  t,  R(t)/R(0)  is  bounded  and  therefore  the  solutions  of  the 
above  equation  are  given  by  (2.9)  (see  Feller  (1968),  p.  459)  □ 

When  A  =  0  in  (2.9),  R(t)  =  R(0)  for  all  t  >  0,  i.e.  equality  holds  in  Holder's 
inequality  (2.6),  and  thus  for  each  pair  s  <  t  we  have  X(t)  =  X(s)  a.s.  Hence  X  is  equal  in 
law  to  a  constant  process  (C(t)  =  aZ,  -oo  <  t  <  oo},  a  >  0,  Z  a  standard  SoS  r.v.,  and  every 
separable  modification  of  X  has  constant  paths.  At  the  other  extreme,  when  A  =  +»,  we 
have  R(t)  =  0  for  t  >  0  and  R(0)  >  0,  so  that  the  stationary  process  X  is  not  continuous  in 
probability  and  thus  its  sample  functions  do  not  have  measurable  modifications  ([2],  p.  3), 
i.e.  it  is  very  irregular.  The  interesting  case  then  is  when  0  <  A  <  oo.  In  the  Gaussian  case 
o=2,  the  symmetry  of  R  and  the  fact  that  it  determines  the  distribution  of  X,  imply  that 
the  only  stationary,  Gaussian,  left  weakly  Markov  processes  are  the  Ornstein-Uhlenbeck 
processes  with  covariance  function  R(t)  =  R(0)e~'^ltl,  -oo  <  t  <  oo,  which  are  in  fact 
Markov.  As  we  shall  see  in  the  non-Gaussian  SoS  case  1  <  a  <  2  there  exist  left  weakly 
Markov  stationary  processes  that  are  not  Markov  (see  e.g.  Section  4). 

Results  analogous  to  Theorem  2.1  and  Corollary  2.1  are  clearly  valid  for  the  right 
weak  Markovian  property.  We  will  not  repeat  the  details  here;  we  only  mention  that  (2.4) 


takes  the  form 


R(tpt2)R(t2,tg)  =  R^pt^R^^)  for  all  tj  <  t2  <  t^, 


and  (2.5)  takes  the  following  form: 


R(t,s)  =  H(s)K(t) 


<a-i> 


for  all  t  <  s, 


where  K(t)/H(t)  is  positive  and  nondecreasing  on  T. 

Therefore,  if  X  is  two-sided  weakly  Markov  with  R(t,t)  i  0  on  T,  then 


(2.11) 


i 


k- 


H1(t)K1(s)<0“l>,  s<t, 

.  H,(t)K2(s)<Q_1>,  t<s. 

When  1  <  a  <  2  the  two  pairs  of  functions  Kp  Hp  and  H2  need  not  be  identical,  as  is 
the  case  with  the  time  changes  of  SaS  Levy  motion  defined  by  (2.7)  where 

(2.12)  H2(t)  =  H1(t)<2-a>K,(t)<a-1>,  K2(t)  =  Hj(t). 

These  Levy  motion  time  changes  are  in  fact  Markov,  as  follows  from 

xf[>  =  +  H(‘)  {Lf^*)) 

and  they  have  o-stable  conditional  distributions  symmetric  about  (H(t)/H(s)}X(s)  for 
s  <  t: 

9  {exp[irX(t)]|X(u),u<s}  =  expfir  X(s)}  Sexp{irH(t)[L(r(t))-L(r(s))]} 

(2.13)  =exp{ir«|i|x(s)-|r|a|H(t)|a[r(t)-r(s)]} 

=  *  (exp[irX(t)]|X(s)}. 

In  particular  every  two-sided  weakly  Markov  stationary  SaS  process  has 

— A,t 

f  R(0)e  ,  t  >  0, 

R(t)  =  \  A9t 

[  R(0)e  ,  t<0, 


for  some  0  <  ApA9  <  00.  When  1  <  a  <  2,  the  exponents  A^  and  Ag  need  not  be  equal,  see 
(2.15). 

Among  the  time  changes  (2.7)  of  SaS  Levy  motion  the  only  stationary  ones  are  of  the 

form 

(2.14)  X(t)  =  ae'At  L(eaAt),  <  t  <  00, 

for  some  0<a<oo,  0  <  A  <  <»;  and  in  fact  it  can  be  easily  seen  they  are  the  only  ones 


J' 


with  stationary  bivariate  distributions  (i.e.  for  these  time  changes,  bivariate  stationarity 
implies  stationarity).  When  0  <  A  <  ®  the  Markov  processes  (2.14)  will  be  called  SoS 
Omstein-Uhlenbeck  with  parameters  a  and  A.  Using  (2.7),  (2.11)  and  (2.12)  we  conclude 
that  the  covariation  function  of  the  SaS  Ornstein-Uhlenbeck  process  (2.14)  is  given  by 


(2.15) 


R(t)  = 


R(0)eAt,  t>0, 
R(0)e^^At,  t  <  0, 


and  is  not  symmetric  unless  a=2. 
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3.  MORE  ON  TIME  CHANGED  LEVY  MOTION 

In  Section  2  we  saw  that  time  changes  of  SaS  Levy  motion  are  Markov  with 
conditional  distribution  of  X(t)  given  X(s)  a-stable  and  symmetric,  for  any  s  <  t.  Here  we 
show  that  all  SaS  Markov  processes  with  right  to  left  conditional  distributions  (as  above) 
o-stable  and  symmetric  are  made  up  from  independent  segments  of  time  changed  Levy 
motion;  and  in  particular  the  only  stationary  ones  with  dependent  values  are  SaS 
Ornstein-Uhlenbeck  processes. 

Recall  that  a  SaS  Levy  motion  L  =  {L(t),  t  >  0}  is  a  process  with  stationary 
independent  increments,  L(0)  =  0  a.s.,  and  for  all  real  r  and  t,  s  >  0, 

(3.1)  8 exp{ir[L(t)-L(s)]}  =  exp{-|r|a  1 1— s | } . 

If  r(t)  is  positive  and  nondecreasing  on  T,  and  H(t)  is  positive  on  T,  then  the  time  change 


of  the  SaS  Levy  motion 


X(t)  =  H(t)  L(r(t)),  t  G  T, 


is  Markov  and  for  s  <  t,  the  conditional  distribution  of  X(t)  given  X(s)  is  a-stable  and 
symmetric,  cf.  (2.13). 

Note  that  we  can  regard  the  above  time  change  as  increasing  (the  new  clock  r(t)  is 
an  increasing  function  on  T).  Similarly,  we  can  define  a  decreasing  time  change  of  SaS 
Levy  motion  by  taking  the  clock  r(t)  in  (3.2)  to  be  a  decreasing  nonnegative  function  on  T. 
Of  course,  the  new  class  of  SaS  processes  obtained  in  this  way  consists  of  Markov 
processes.  Moreover,  they  have  the  following  common  property  for  s  >  t,  the  conditional 
distribution  of  X(t)  given  X(s)  is  o-stable  and  symmetric.  These  properties  of  the  time 
changes  of  SaS  Levy  motion  are  quite  remarkable.  We  will  see  in  this  section  that  there 
are  not  many  Markov  SaS  processes  whose  conditional  distributions  are  a-stable  and 
symmetric. 
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The  conditional  distributions  of  every  Gaussian  process  are  Gaussian  and  symmetric 

(around  the  conditional  mean).  Non-Gaussian  stable  processes  in  general  do  not  have 

stable  conditional  distributions.  Our  aim  in  this  section  is  to  characterize  the  classes  ^ 

a 

and  of  all  SaS  Markov  processes  X  which  have  the  property  that  for  all  s  <  t  (s  >  t, 
correspondingly)  the  conditional  distribution  of  X(t)  given  X(s)  is  o-stabie  and  symmetric. 
Of  course  if  1  <  a  <  2  and  ,/{X(t)  |  X(s)}  is  symmetric  about  some  point,  this  point  of 
symmetry  is  necessarily  the  conditional  mean  £{X(t)  |X(s)}. 

Theorem  3.1  characterizes  the  processes  in  and  when  1  <  a  <  2  Those 
with  covariation  function  nonvanishing  everywhere  are  time  changes  of  Levy  motion,  as  in 
(3.2).  The  general  process  in  (jf!^)  is  then  made  up  of  independent  segments  of 
Levy  motion  time  changes  on  disjoint  intervals.  There  are  two  extreme  (and  uninteresting) 
cases:  the  (very  smooth)  constant  process  (X(t)  =  Z  a.s.  for  each  t),  and  the  (very  rough) 
process  consisting  of  independent  random  variables.  In  the  stationary  case  the  latter 
process  would  have  independent  and  identically  distributed  SaS  r.v.'s,  with  scale  parameter 
a  >  0,  and  we  denote  it  by  I  =  (I  (t),  -  ®  <  t  <  »}.  The  only  stationary  processes  in 

a  a. 

are  the  SaS  Ornstein-Uhlenbeck  processes  (2.14)  and  those  with  iid  values:  I&.  The 
only  stationary  processes  in  are  the  inverted  SaS  Ornstein-Uhlenbeck  processes 


defined  by 


X(t)  =  aeAt  L(e  aAt),  -«  <  t  <  ®,  a  >  0,  A  >  0, 


and  the  processes  I  .  The  SaS  Ornstein-Uhlenbeck  processes  (2.14)  and  (3.3)  coincide 

a 

trivially  in  the  Gaussian  case  o=2,  but  not  in  the  case  a  <  2  (more  on  this  point  is  said  in 
Theorem  3.2). 

In  the  statement  of  Theorem  3.1  equality  in  law,  ,  means  equality  of  all  finite 
dimensional  distributions. 
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Theorem  3.1.  Let  X  belong  to  ( correspondingly ,  JC^)  for  some  1  <  a  <  2. 

(i)  //  its  covariation  function  satisfies  R(t,s)  i  0  for  alls  <  t  in  T,  then 


{X(t),t€T}  =  (H(t)L(r(t)),  t  €  T} 


for  some  positive,  nondecreasing  ( correspondingly ,  nonincreasing)  function  r  on  T,  and 
some  positive  function  H  on  T. 

(it)  If  X  is  stationary  then  either 


(X(t),  —oo  <  t  <  oo}  =  {ae  ^L(ea/^),  —  oo  <  t  <  oo} 


for  some  a  >  0  and  0  <  A  <  oo  ( correspondingly ,  0  <  -A  <  oo),  or  else  for  some  a  >  0, 


{X(t)  ,  -00  <  t  <  X 

CL 

In  order  to  prove  Theorem  3.1  we  need  the  following  properties  of  bivariate  SaS 
distributions  which  are  of  independent  interest.  Let  us  recall  that  the  r.v.'s  X^  and  X2  are 
jointly  SaS  if  their  joint  characteristics  function  is  of  the  form 

frexpfiOrjXj+^Xj)}  =  exp{  -fs  |r1x1+r2x2|ftdr(x1,x2)  } 

2 

for  all  real  ,r2,  where  F  is  a  uniquely  determined  symmetric  finite  measure  on  the  unit 
o 

circle  S2  in  IR  .  When  1  <  a  <  2,  the  covariation  of  X^  with  X2  is  given  by 

Cov[Xp  X2]  =  /g  x1x2<Q^1>dr(x1,x2) 

[5]  (which  is  consistent  with  (2.2)).  We  denote  by  ||Xj||  their  scale  parameter  ||Xjj|^  = 
/ g  |  Xj  |  adr(Xj,x9)  =  Cov[Xj,  X-] ,  i=l,2,  and  we  have  by  Ranter  (1972), 


X 


^(X2|X1)  = 


Cov[X2,  Xj]  A 


xjt  xi  ~  p2ixr 


Proposition  3.1.  Let  X^  and  X2  be  jointly  SoS  and  1  <  a  <  2.  Then  the  following  are 
equivalent. 


(i)  J^(X9|Xj)  is  ot-stable  and  symmetric. 

( ii )  _^2l^l 's  independent  o/Xy 

(in)  f  is  concentrated  on  ±  (0,1),  ±  ((1+/>21)-1/2,  />21(l+p21)_1//2). 


Under  any  of  these  conditions  we  have 


(3.4)  IIX^XX  =  l|X2||“-  I/.21|“||X1||“  =  |X2I|"  -  . 


Xj,X9  are  independent  iff  Cov[X2,  Xj]  =  0  iff  /?21  =  0. 


Proof.  Assume  (i).  Then 


*  {exp(ir2X2)|Xj}  =  exp{-|r2|aM(X1)  +  ir2N(Xj)}, 


for  some  real  measurable  functions  M  and  N  with  M  >  0.  It  follows  that 


N(X1)=  »(X2|X1)  =  p2lXr  Let 


Z  =  p,1X1  +  M(X1)1''aZ0 

SJu  1 

where  Zq  is  independent  of  Xj  and  8  exp(irZQ)  =  exp(-|r|  ),  r  6  K  .  Then  clearly 


./(Z|Xj)  =  J<(X2|X^)  and  thus  Jf(XyZ)  =  «/(Xj,X2).  It  follows  that  Z  -  /^X^  is 
SaS,  and  thus  for  some  c  >  0  and  every  real  r, 


exp(-|r| ac)  =  JfexpfirfZ-^jXj)}  =  Sexp  {irM(X1)1/°Z()}  =  8 exp{-|r| "M^)}. 


By  the  uniqueness  of  the  Laplace  transform  we  conclude  that  M(Xj)  =  c  a.s..  Then 


.v  *•  -  V  vvvvvvv;/..-.  «%  ‘ 


Z“^21X1  =  c^0  'ls  '^e^ndent.  of  Xp  and  this  implies  that  ^2-/?21Xl  *s  a^so 
independent  of  Xp  proving  (ii). 

Conversely,  assuming  (ii)  and  writing  X2  =  (X^p^Xj)  +  P2iXi  we  obtain’  since 
X2  p2i Xj  isS aS, 

(3.5)  $  (exp(ir2X2)  f  }  =  exp{—  (r2(  ll^2_^21Xl^a  *r2^21Xl^ 

Ck 

so  that  (i)  is  satisfied,  an^  in  fact  the  constant  M(X^)  =  c  is  equal  to  ll^2~^21Xl^a 
also  obtain 

I|X2||“  =  ||X2-/.21X1||“+  |/)21|a||X1||“ 

from  which  (3.4)  follows. 

Therefore  (i)  or  (ii)  imply 

(TexplifrjXj+rjXj))  =  JTexpfi^+r^pX,  -  |r2|Q||X2-/>21X1||“) 

=  exp H r j  +r2P2 j I "||X , || “ -  | r2 1“|| Xj-p2 , X, || "} 
from  which  (iii)  is  evident.  And  conversely,  assuming  (iii)  we  have 

^exp{i(r1X1+r2X2)}  =  exp{-|r2|  “dj  -  |rj+r2P21 1  %} 
for  some  dpd2  >  0,  and  therefore 

(S  (i[a(X2-p21X1)  +  bXj}  =  8  exp{i[(b-ap21)X1  +  aX>2]} 

=  exp{-|a|ad1  -  |b|ad2} 

from  which  (ii)  follows.  D 

Corollary  3.1  Let  X1  and  X9  be  jointly  SaS  and  1  <  a  <  2.  Then  the  following  are 
equivalent. 


(i)  J/(X2|X^)  and  Jr(X^|X0)  are  ot-stable  and  symmetric. 

(it)  X1  and  are  either  independent  or  linearly  dependent. 

(Hi)  f  is  concentrated  on  {±(0,1),  ±(1,0)}  or  on  ±  (c,  (1-0^)^“)  for  some 


0  <  c  <  1. 

Proof.  In  view  of  Proposition  3.1,  (i)  implies  that  T  is  concentrated  on  the  set  {±(0.1), 

±((l+/>2 1)-1^  f,21^1+p21^_1^2^  and  also  on  the  set  W1’0)’  ±(p19(l±^2)_1//2, 
(l±p29)-1^)},  from  which  (iii)  follows.  The  converse  is  clear,  as  is  the  equivalence  of  (ii) 

and  (iii).  □ 

Proof  of  Theorem  3.1 .  Because  of  symmetry  we  will  consider  only. 

(i)  Since  for  all  s  <  t  in  T,  J^{X(t)  |  X(s)}  is  ce-stable  and  symmetric,  it  follows  from 
Proposition  3.1  (cf.  (3.5))  that 

(3.6)  *{exp[irX(t)]|X(s)}  =exp{-|r|a||  X(t)-ptgX(s)||“  +  irptsX(s)} 

where  ^  =  R(t,s)/R(s,s)  and  ||X(t)-^tsX(s)||  J  =  R(t,t)  -  | R(t,s) |  “/R^s)*"1  (cf.  (3.4)). 
Since  X  is  Markov  and  R(t,s)  f  0  for  s  <  t,  by  Theorem  2.1,  R(t,s)  has  the  representation 
(2.5),  which  implies 


^t,s  H(s)  ’ 

||X(t)-(>tsX(s)||“=  H(t)K(l)<0“1>-  |H(t)|"(Sf|)°“1  =  |H(t)|a{r(t)-r(s)}. 


where  r(t)  =  {K(t)/H(t)}ft_  .  Thus 


H(t 


8  {exp[irX(t)]|X(s)}  =  exp{-|r|  |H(t)|  [r(t)-r(s)j  +  ir  H(s)  X(s)} ’ 

which  is  the  same  as  the  conditional  characteristic  function  for  Y(t)  =  H(t)L(r(t)),  t  e  T. 
given  in  (2.13).  It  follows  that  y{X(t)|X(s)}  =  Jf{Y(t) | Y(s)}  for  all  s  <  t.  It  is  also 
clear  from  (2.S)  that  ||X(t)||®  =  ||Y(t)H*,  so  that  ^{X(t)}  =  y{Y(t)}  for  all  t,  and  since 


K< 


SWWi'3 


S  V 


both  X  and  Y  are  Markov,  Jf(X)  =  Jf(Y).  The  function  H  may  be  taken  positive  without 
loss  of  generality. 

(ii)  If  X  is  stationary,  by  Corollary  2.1,  R(t,s)  =  R(0)e”'^t_s)  for  all  s  <  t,  for  some 
0  <  A  <  oo.  When  A  =  +  oo,  R(t,s)  =  0  for  all  s  <  t,  and  by  Corollary  3.1  it  follows  that 
./(X)  =  for  some  a  >  0.  When  0  <  A  <  ®,  part  (i)  of  the  theorem  is  applicable, 

necessarily  with  H(t)  =  ae-/^  and  r(t)  =  bea/^  for  some  a,b  >  0.  This  completes  the 
proof.  o 

It  is  worth  mentioning  that  in  the  Gaussian  case  (o=2),  clearly  JtL^  =  In 

contrast  when  1  <  a  <  2  the  common  processes  of  and  are  few  and  trivial. 

O'  a 

Theorem  3.2  (i)  There  is  a  one-to-one  coirespondence  between  and  given  by 

X(t)  =  Y ( r ( t ) ) ,  t  G  T,  X  6  ^0^ •>  ^  e  ’  for  any  f uncti°n  T:  T  ->  T  one-to-one , 
onto ,  and  such  that  r(s)  >  r(t)  if  s  <  t  (e.g.  r(t)  =  -t  when  T  =  R*). 

(ii)  A  process  X  belongs  to  n  1  <  a  <  2,  if  and  only  if  it  is  of  the 

following  form .  for  some  finite  or  denumerable  set  of  intervals  in  T,  X(t)  = 

a(t)Xn  a.s.  for  each  t  €  I  ,  n=l,...,N,  where  a  is  a  real  function,  nonvanishing  in  the 
interior  of  and  the  r.v.'s  {(Xn)^_j,  X(t),  t  6  T  \  1^}  are  independent. 

(Hi)  The  only  stationary  processes  in  JC^  fl  1  <  a  <  2,  are  the  processes  I 

with  iid  r.v.'s  and  the  constant  processes  (X(t)  =  aZ,  <  t  <  oo,  a  >  0,  Z:  SaS  r.v.). 

Proof,  (i)  is  clear. 

(ii)  Suppose  X  €  n  Then  by  Corollary  3.1,  for  each  s,t  in  T,  the  r.v.'s 

X(s)  and  X(t)  are  either  independent  or  linearly  dependent.  But  since  X  is  Markov,  if  X(t) 
is  a  multiple  of  X(s),  it  will  necessarily  be  a  multiple  of  each  X(u)  for  u  in  between  s  and  t. 
Hence  for  each  t  6  T  there  is  an  interval  1^  such  that  for  all  s  6  1^,  X(s)  is  a  multiple  of 


X(t).  Denoting  by  {In}^[_i  those  intervals  among  the  It,  t  6  T,  with  positive  Lebesgue 
measure  we  obtain  the  result. 

(iii)  In  view  of  stationarity,  if  two  r.v.'s  of  the  process  are  linearly  dependent,  they 
will  all  be  a.s.  equal;  and  if  two  r.v.'s  of  the  process  are  independent,  they  will  all  be 
independent.  D 
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4.  SUB— GAUSSIAN  PROCESSES 


Another  important  class  of  SaS  processes  is  the  class  of  sub-Gaussian 
processes  which  is  defined  by 

(4.1)  X(t)  =  A1/2G(t),  teT, 

where  G  =  {G  ( t ) ,  teT}  is  any  Gaussian  process  and  A  is  a  totally  skewed  to  the 
right  (a/2)-stable  random  variable  (see  Feller  (1968))  which  is  positive  with 
probability  1  and  satisfies 

(4.2)  £ exp{-uA}  =  exp{-ua/2} ,  u  >  0. 

A  sub-Gaussian  process  (4.1)  is  easily  seen  to  be  a  SaS  process,  and  when  T  =  IR*, 
X(t)  is  stationary  if  and  only  if  G(t)  is. 

Our  task  in  this  section  is  to  investigate  if  we  can  find  a  Markov,  or  at  least 
weakly  Markov  process  among  sub-Gaussian  processes.  The  obvious  candidates  are 
the  sub-Gaussian  processes  given  by  (4.1)  with  G  being  Gaussian  Markov  process. 
For  example,  G(t)  could  be  a  Brownian  motion;  in  this  case  we  call  the 
corresponding  process  X(t)  a  sub-Brownian  motion.  If  T  =  IR^  and  G(t)  is  Ornstein- 
Uhlenbeck  process,  we  call  X(t)  sub-Omstein-Uhlenbeck  process.  The  first  result 
characterizes  the  weakly  Markov  sub-Gaussian  processes. 

Theorem  4.1.  Let  X  be  a  sub-Gaussian  SaS  process  given  by  (4.1).  Then  the  following  are 
equivalent. 

( i)  X  is  left  weakly  Markov. 

( ii )  X  is  right  weakly  Markov. 

( iii )  G  is  Markov. 


Proof.  We  note  for  future  reference  the  following  fact  (see  [5] ) :  If  (Yj  ,Y2)  is  a 
zero-mean  Gaussian  vector  in  R2  with  <$  Y2  =  $  Y2  =  1 ,  S  Y^Y 2=  p,  and  if  X^  =  2Yj , 
X2  =  A1/^Y2,  where  A  is  distributed  according  to  (4.2)  and  is  independent  of 
(YpY2),  then 

(4.3)  Cov[Xj,  X2]  =2"®/V 


Assume  (i) .  Then  by  Theorem  2.1  the  covariation  function  of  X  satisfies  (2.4) 
and  by(4.3)  we  conclude  that  the  covariance  function  of  the  Gaussian  process  G  in 
(4.1)  satisfies  the  same  relation.  Thus  G  must  be  Markov  (since  it  is  Gaussian), 
i.e.  (iii)  holds. 

Conversely,  assume  (iii);  we  will  show  ( i) .  For  any  s  <  t ,  any  u^ ,u2, . • .  ,u.  <  s 
and  any  a^  ,a2, . . .  ,a,  ,  we  have  by  (4.3) , 


Covtx<t>-|fel}x(s>-Ei=iaix(ui)] 

=  Cov[A^{C(t)  -MM5M)  C(s)},  A!/24 , 4,0(0,)] 

S(G^(s))  1-1  1  1 

=  2"a/2[^_1  aiG(ui)}2](Q^2)/2  Cov [G (t )  G(s),  a.G(u.)]  =0 

£(G  (s)) 


since  G  is  Markov.  Therefore,  by  Theorem  2.1,  X  is  left  weakly  Markov,  i.e.  (i) 
holds.  The  equivalence  of  (ii)  and  (iii)  is  now  deduced  by  the  symmetr>  .  □ 


Corollary  4.1.  Any  sub-Gaussian  weakly  Markov  SaS  process  with  covariation  function 
R(t,s)  t  0  for  any  t,s,  is  equivalent  to  a  nondecreasing  time  change  of  sub~Brownian 
motion. 


Proof .  Formula  (4.3)  implies  that  the  covariation  function  of  a  sub-Gaussian 
process  determines  its  finite  dimensional  distributions.  Therefore,  we  have  only 


to  demonstrate  that  the  covariation  function  of  a  sub-Gaussian  weakly  Markov  SaS 
process  can  be  realized  by  a  time  changed  sub-Brovnian  motion.  This  is  easy.  By 
Theorem  2.1  and  (4.3),  the  covariation  function  of  a  sub-Gaussian  weakly  Markov 
SaS  process  can  be  represented  as  R(t,s)  =  H(t)K(s)<Q‘_*>  for  all  s,t,  where 
K(t)/H(t)  is  positive  and  nondecreasing  on  T.  Let  G^(t)  =  H(t)B(r(t)),  where  B  is 
the  standard  Brownian  motion,  and  r(t)  =  2[K(t)/H(t)]^Qf-^/<*.  It  is 
straightforward  to  check  using  (4.3)  that  the  process  X^(t)  =  A^G^(t)  has  R(t,s) 
as  its  covariation  function.  Since  is  a  time  change  of  sub-Brownian  motion,  the 
proof  is  complete.  □ 

Among  all  SaS  sub-Gaussian  processes  only  very  few  and  highly  degenerate  are 
Markov.  All  these  processes  are  characterized  in  the  following  simple  fashion.  We 
say  that  a  deterministic  function  a(t),  t  e  T,  is  bom  at  the  level  a^  and  dies  at  the 
level  if  for  some  <  S2 »  which  may  be  boundary  points  of  T,  a(t)  =  a^  for  t  <  s^ 
(or  t  <  Sj) ,  a(t)  =  a2  for  t  >  s2  (or  t  >  s2) ,  and  a(t)  i  {a1?a2}  outside  of  the  above 
intervals;  s^  is  the  birth  time  and  s2  is  the  death  time. 

Theorem  4.2.  A  SaS  sub-Gaussian  process  given  by  (4.1)  is  Markov  if  and  only  if  the 

Gaussian  process  G  has  one  of  the  following  forms :  G(t)  =  a(t)  or  G(t)  = 

a(t){Y1  l(t  <  tQ)  +  Y2  l(t  >  tQ)}  or  G(t)  =  a(t){Y1  l(t  <  tQ)  +  Y2  l(t  >  tQ)},  t  e  T, 

2  2  2 

where  (Y^,Y2)  is  a  jointly  Gaussian  vector  in  !R  ,  tfY^  =  £Y2  =  1,  tg  e  T,  and  a(t)  is  a 
real  function  that  is  bom  and  dies  at  the  level  zero. 

A  technical  result  precedes  the  proof  of  the  theorem. 

Lemma  4.1.  Let  Gj,  G2,  G^  be  i.i.d.  standard  normal  random  variables  and  let  W  be  a 
positive  random  variable  not  equal  a.s.  to  a  constant  and  independent  of  G^,  G2,  G^. 
Then  there  is  a  Borel  set  B  such  that  P(WGj  €  B)  >  0  and  for  any  Xj  6  B  the  random 
variables  WG2  and  WG^  are  not  independent  given  WG^  =  x^. 


24 


Proof.  Suppose  on  the  contrary  that  for  almost  every  value  of  WG^  the  r.v. ' s 
WG^  and  WG.-,  are  independent  given  WGj  =  x^.  Then  for  almost  any  Xp  x2,  we  have 


(4.4) 


f i(y !xi)  =  f2(ylxi’x2) 


for  almost  any  y,  where  fj(-  |x^)  is  the  conditional  density  of  WG3  given  WG1  =  Xj, 
and  f2(-  IxpXg)  is  the  conditional  density  of  WG3  given  WGj  =  x2,  WG2  =  x2>  We  have 

/o(2?r)“1/2w_2exp[-(y2+x^)/(2w2)]dF(w) 


(4.5) 


f ! (y |xj)  = 


(4.6)  f2(y|x1,x2)  = 


/o  w  exp[-x^/(2w2)]dF(w) 

/J(27r)_1//2w~3exp[-(y2+xJ+x2)/(2w2)]dF(w) 
!q  w~2exp[-(x2+x2)/(2w2)]dF(w) 


where  F  is  the  distribution  function  of  W.  It  follows  from  (4.5)  and  (4.6)  that 
f j  ( - )  * )  is  continuous  on  K  ,  and  f2( • | • , • )  is  continuous  on  Or.  We  conclude  that 

(4.4)  is  equivalent  to 

(4-7)  g3(y2+x2+x2)  g^x2)  =g2(y2+xi)  g2(xi+x2)’ 

where 

gn(r)  =  /q  w~nexp[-r/(2w2)]dF(w),  r  >  0. 

By  Holder's  inequality  we  conclude  that  the  left  hand  side  of  (4.7)  is  strictly 
larger  than  its  right  hand  side  for  all  triples  (y,Xj  ,x2)  of  the  kind  (0,Xj  ,0) .  By 
the  continuity  we  conclude  that  there  is  an  e  >  0  such  that  (4.7)  and,  therefore, 

(4.4) ,  do  not  hold  for  the  triples  (y,x^,x2)  €  (-e,t)  .  This  contradicts  the 

assumption  that  for  almost  any  x^  and  x2,  (4.4)  is  true  for  almost  any  y.  a 

Proof  of  Theorem  4.2.  Suppose  that  the  sub-Gaussian  process  X  given  by  (4.1)  is 
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Markov.  If  all  r.v.'s  of  the  Gaussian  process  G  are  linearly  dependent,  then 
G(t)  =  a(t)Yj  for  some  standard  normal  r.v.  Y^.  Now  assume  G  has  two  linearly 
independent  r.v.'s,  say  G(tj)  and  G(t2),  tj  <  t2-  Furthermore,  assume 
£{G(tj)G(t2)}  ^  0,  as  the  argument  is  even  simpler  in  the  case  of  independence. 
Then  we  can  write  G(t,2)  =  a^Gj,  G(tj)  =  a^Gj  +  a12G2,  where  Gj,  G2  are  i.i.d. 
standard  normal  r.v.  's  and  the  coefficients  a^,  aj2»  a2l  are  different  from  zero. 
Take  any  t^  >  t2.  Then  we  can  write  G(t^)  =  a^G^  +  ag2G2  +  a^G^,  where  G^  is  a 
standard  normal  r.v.  independent  of  (Gj,G2).  The  process  G  is  Markov  by  Theorem 
4.1;  therefore  its  covariance  function  satisfies  the  relation  (2.4).  Rewriting 
this  relation  for  the  triple  (t^t^tg)  in  terms  of  our  particular  representation 
ofG(t-),  i= 1 ,2,3,  we  obtain 

(alla31  +  a12a32^  a21  =  (ai ia2 1 ) (a21a3l) ’ 

and  since  a^,  a^2  and  a21  are  different  from  zero,  we  conclude  that  a^2  =  0. 

Recall  now  that  by  the  assumed  Markovianness  of  X(t) ,  for  almost  any  x2,  X(t^) 
and  X(tg)  are  independent  given  X(t2)  =  x2-  Since  a^2  =  0,  it  follows  that 
aj2A^G2  and  a-^W^Gg  are  independent  given  a^A^Gj  =  x^.  Suppose  first  that 
a^g  t  0.  Then  it  follows  that  for  almost  any  Xj,  A^G2  and  A^Gg  are  independent 
given  A^^Gj  =  Xj/a2j .  But  this  is  impossible  by  Lemma  4. 1 .  Therefore  agg  =  0.  It 
follows  that  G(tg  )  -  agjGj,  and  thus  for  any  t  >  t2,  G(t)  is  linearly  dependent  on 
G(t2).  The  same  argument  shows  that  for  any  t  <  t^,  G(t)  is  linearly  dependent  on 
G(tj),  and  that  there  is  a  tQ  €  [tj,t2]  such  that  for  any  t^<  t  <  t  ,  G(t)  is 
linearly  dependent  on  G(tj),  while  for  any  t  <  t  <  t2,  G(t)  is  linearly  dependent 
on  G(t>2).  G(tQ)  itself  has  to  be  linearly  dependent  either  on  G(t^)  or  on  G ( 1 2 )  ■ 
This  proves  the  dependence  structure  of  G(t)  described  in  the  statement  of  the 
theorem. 


kj  uppvjov.  nuw  tuau  ^  pv/xxtoo  o  ^  x  x  ouou  iiiiat  a.  ^  k?  ^  /  T  C*''  2'  ~  '"/ 

and  a(s^)  ^  0.  Applying  the  assumed  Markovianness  of  X(t)  to  the  triple  (s^^jSg) 
shows  that  one  of  the  following  pairs  of  r.v.  's  is  independent:  (A^2Yp  A*/2Y^)  or 
{KXlhv  KX!2\2)  or  (A1^2Yj  ,  A^2Y2).  The  r.v.’s  of  the  first  two  pairs  are 
obviously  dependent,  and  those  of  the  third  pair  are  also  known  to  be  dependent 
(see  Lemma  2.1  in  [6]).  Therefore  a(t)  must  be  born  and  die  at  the  level  zero. 
Since  it  is  obvious  that  for  any  Gaussian  process  of  the  above  form  the 
corresponding  sub-Gaussian  process  is  Markov,  the  proof  of  the  theorem  is 
complete.  o 

In  the  particular  case  of  a  stationary  sub-Gaussian  SaS  process  we  can  draw 
the  following  simple  conclusion  from  the  above  results. 

Corollary  4.2  The  only  weakly  Markov  stationary  sub-Gaussian  SaS  processes  are  the 
sub-Omstein-Uhlenbeck  processes  and  the  constant  processes.  The  only  Markov  stationary 
sub-Gaussian  SaS  processes  are  the  constant  processes. 

Corollary  4.2  shows  that  the  sub-Ornstein-Uhlenbeck  process  is  the  only 
weakly  Markov  stationary  sub-Gaussian  process.  However,  the  class  of  sub-Gaussian 
processes  is  not  closed  under  linear  combinations  of  its  independent  members. 
Therefore  a  natural  question  arises:  could  we  obtain  a  new  weakly  Markov 
stationary  SaS  process  as  a  linear  combination  of  independent  stationary  sub- 
Gaussian  processes? 

Let  Gj,  G2 * • • - *G  be  independent  stationary  Gaussian  processes  such  that  for 
any  t, 

(4.8)  ^{Gi(t)Gi(0)}=/?i(t),  ^(0)  =  1,  i=l ,2, . . .  ,n. 

We  assume  that  all  the  correlation  functions  p^,p2, • ■ • ,Pn  are  different.  Let 


A-, , A0 , . . . , A  be  i.i.d.  random  variables  distributed  according  to  (4.2)  and 

i  z  n 

independent  of  the  Gaussian  processes  Gj,G2, . .  ■  ,G  .  Let  b^ , b2 » •  •  •  ,b0  be  positive 
real  numbers.  We  will  prove  that  if  n  >  1,  the  stationary  SaS  process 


(4.9) 


X(t)  *£j=1  biAV2Gi(t),  -®  <  t  <  ao, 


cannot  be  weakly  Markov.  We  start  with  the  following  lemma. 

Lemma  4.2.  Let  S  ->  R* ,  i=l,2,...,n,  be  arbitrary  distinct  functions  on  a  set  S. 
Then  there  is  a  finite  subset  {Sj^,  •  •  •  ,sk)  C  S  (k  <  n-1)  and  real  numbers  ^ , 02, . . . , 0j. 
such  that 

Ej_l  ^ j ^1  ( s j )  *  Sj=l  ^jV»i(sj) *  any  i=2’--->n- 

Proof.  We  prove  the  lemma  by  induction  in  n.  For  n=2  the  claim  of  the  lemma  is 
trivial,  and  one  can  choose  k=l  and  =  1.  Suppose  that  the  claim  of  the  lemma  is 
correct  for  n  =  nQ.  We  shall  prove  this  claim  for  nQ+l  functions,  •  •  • , <Pn  , 

ipa  +j.  By  the  assumption  of  the  induction  we  know  that  there  is  a  set 

{Sj ,s2,  •  • .  ,sk  }  C  S,  and  real  numbers  0j,02>  •  •  •  »  such 

k  k 

A(i)  =  ^(Sj)  ^j^j(Sj)  t  0,  any  i=2, . . .  ,nQ. 

If  A(nQ+l)  t  0,  then  there  is  nothing  to  prove,  so  we  assume  that  A(nQ+l)  =  0. 

Then,  there  is  an  sk  +J  e  S  such  that  ^(s^  +1)  i  <pn  +1(sk  +1).  Put 

7i  =  ^i^skQ+l^  -  (skQ+l) ’  1=2’---’n0- 

Clearly,  we  can  choose  a  real  number  #k  +1  satisfying  the  following  conditions: 

Ko 


(b)  0^  +j  ^  -A^/7^  for  all  such  i=2, ...  ,nQ,  for  which  7.^  0. 


k  +1  k  +1 

Sj°i  (sj)  #  sj=i  ^j^i(sj)  for  any  i=2,...,nQ,  nQ+l. 

so  that  the  proof  of  the  lemma  is  complete.  □ 

Proof  of  the  claim.  If  the  process  given  by  (4.9)  were  weakly  Markov,  we  would 
have  for  any  t ,  any  r  >  0 ,  any  s  ^  >  0 ,  i=l , . . . , k ,  any  real  numbers  c ,  9^ , . . . ,  A  , 
that  for  some  0  <  A  <  00, 

(4.10)  Cov[X(t+r)  -e~ArX(t),  X(t)  0-X(t-s-)]  =0. 

J  J  J 

Since  the  sub-Gaussian  processes  A^G^t),  i=l,...,n,  are  independent,  we 
conclude  by  the  properties  of  covariation  (see  Veron  (1984))  and  by  (4.3)  and 

(4.10)  that 

0  =  e;=1  bf  Cov[Al/2{G.(t+r)  -  ApfOjd)  ♦  c  ^(t-sp}] 

=  2 -°/2  S"=1  bf  f Var{Gj (t )  +  c  5*=1  ^Cjft-Sj)}]  (q“2)/2  * 

(4.11)  Cov[Gi(Ur)-e-ArGi(t),Ci(t)+c^,  tfjCi(t-s_})] 

=  2~a/2  S?=1  bf  (1  *  2c  »j(..(5j)  *  c2  Ej=1  4,  Yf^Sj-Sf)]  W/K 
[<?;(  r)-*TAr)  *  c  Ej=1  -  eAr  ^|(Sj)}] . 

That  is,  for  any  t ,  any  r  >  0,  any  s^  >  0,  i=l , . . .  ,k,  any  real  c,  0^ , . . .  ,^,  we  have 


m 


■v'« 
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cSi=l>l  »jf/>i(Sj+r)-e'A,'oi(sj)}bi 

(4.12)  x  [1  +  2c  Ej_j  #/i(sj)  *c2Sj=1E^=1  tfip i(sj-sf)]<°^2>/2 

=  5?=1  b?  {e_Ar-/>i(r)}[l  +  2c  5*=1  ^(sp  *  c2  sj=1  £^,  (°“2)/2. 

Since  the  correlation  functions  p^,...,p  are  all  different,  we  can  find  a  t 
such  that  not  all  the  expressions  e~^T-p^  (r) ,  i=l, . . .  ,n,  are  equal  to  zero.  We  can 
assume,  without  loss  of  generality,  that  e  -p^(r)  f  0  for  i=l,...,n  ,  and 
e  =  0  for  i=nQ+l, ...  ,n,  for  some  1  <  nQ  <  n.  According  to  Lemma  4.2,  there 

are  points  s ^  >  0 ,  and  real  numbers  8- ,  i=l , . . . , k ,  such  that 

Sj=i  Vl(Sj)  ^  ^=1  ^/i(sj)’  i=l , . . .  ,n . 

Then  for  these  fixed  values  of  r,Sj, . . .  ,s^,  we  can  rewrite  (4.12)  as 

(4.13)  cE^1a(1)(l  +  2c/ii  +  c2a?)(^/2  =  e“°j  a(2) (1  +  2cMi+c2a?) (^2)/2 

for  certain  real  numbers  a>^  ,  p^,  cr  ,  i=l, . . .  ,n,  a|2^ ,  i=l, . . .  ,nQ.  Then  the  fact 
that  for  any  0  <  d  <  1/2,  any  p,  the  following  family  of  functions  of  c, 


{  (1  +  2c//i  +  c2^)^,  c(l  +  2c/r  +  cVj)-*,  i=l,...,p} 

o 

is  linearly  independent  as  long  as  all  the  pairs  (/z-,^)  are  different,  implies 
that  a|2^  =  0.  But  we  know  that  aj2^  f  0.  This  contradiction  proves  that  no  process 
of  the  kind  (4.9)  can  be  weakly  Markov.  □ 


2  2,-d  : 


4« 
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5.  MOVING  AVERAGES 


A  SaS  moving  average  is  a  process  of  the  form 


X (t )  =  /°°  f  (t— s)  dL (s) ,  —oo  <  t  <  oo , 


where  L  is  an  independeeeently  scattered  SaS  measure  on  (IR^ ,  ^  ,Leb) ,  i.e.  {L(s), 
-oo  <  s  <  oo}  is  a  SaS  Levy  motion,  and  f  e  La(K^ ,  3^  ,Leb) .  X(t)  is  clearly  a 
stationary  process.  When  f  vanishes  on  the  negative  line,  X  is  a  nonanticipating 


moving  average  and 


(5.1) 


X(t)  =  /\f(t-s)dL(s). 


A  nonanticipating  moving  average  is  called  invertible  if  sp{X(t),  t<r}  = 
A  _ 

sp{AL(t) ,  t<r}  =  sp{L(t)-L(s) ,  s<t<r}  for  all  r,  i.e.  the  increments  of  L 
represent  the  innovations  of  X. 

We  now  determine  all  nonanticipating,  invertible  SaS  moving  averages  which 
are  left  weakly  Markov;  these  are  in  fact  Markov,  and  we  shall  show  that  they  are 
precisely  the  SaS  Ornstein-Uhlenbeck  processes,  thus  obtaining  an  integral 
representation  for  these  processes. 

Theorem  5.1  X  is  a  left  weakly  Markov ,  nonanticipating ,  invertible,  SaS  moving  average 
with  1  <  a  <  2 ,  if  and  only  if  it  is  of  the  form 


(5-2) 


X(t)  =  a  /l  e  s)dL(s),  ^>  <  t  <  ®, 


for  some  0  <  a,  A  <  oo,  if  and  only  if  it  is  a  SaS  Ornstein-Uhlenbeck  process. 

Proof .  Let  X  be  a  nonanticipating,  invertible  SaS  moving  average  as  given  by 
(5.1).  Then  X  is  a  left  weakly  Markov  iff  for  some  0  <  A  <  oo, 

Cov[X(r)  -e~ArX(0),  Y]  =  0 


for  all  r>0  and  all  Y  e  sp{X(t) ,  t<0}  =  sp  {AL(t) ,  t<0} ,  i.e.  forallY  =  /^gdL 
where /^Jg(s)  |  ads  <  ®;  i.e.  iffVr>0,  V  g  e  LQ;( (— oo,0) ,  <#(-oo,0)  ,Leb) , 

0  =  Cov[/^f(r-s)dL(s)-e-Ar/^f(-s)dL(s),  J^g(s)dl(s)] 

-  f(r-s)g(s)<Q“1>ds  -  e-Ar/^f(-s)g(s)<0“l>ds. 

Now  putting  g(s)  =  l^_x  (s) ,  x  >  0,  we  obtain  that  if  X  is  left  weakly  Markov  then 
for  all  r  >  0,  x  >  0, 

/®xf(r-s)ds  =  e'Ar/“xf(-s)ds, 

i.e. 

/"xf(u)du  =  e_Ar  /q  f  (u)du, 
and  putting  F(x)  =  / g  f(u)du, 

F(r+x)  =  F ( t)  +  e~ArF(x),  V  r,  x  >  0. 

The  parameter  A  cannot  take  the  values  0  or  +®.  Indeed,  if  A  =  +oo  then 
F(r+x)  =  F(r) ,  V  r,  x  >  0,  implies  F(x)  =  Const. ,  x  >  0,  hence  f  (x)  =  0  a.e.  on  (0,®) 
and  X(t)  =  0  a.s.  for  all  t,  i.e.  the  process  X  is  identically  zero.  On  the  other 
hand,  if  A  =  0  then  F(r+x)  =  F(r)  +  F(x) ,  V  r,  x  >  0,  implies  F(x)  =  cx,  x  >  0,  hence 
f(x)  =  c  a.s.  on  (0,®)  which  does  not  belong  to  Lft(Leb)  unless  c  =  0  in  which  case 
again  the  process  X  is  identically  zero.  It  follows  that  0  <  A  <  ®.  Interchanging  t 
and  x  we  also  have 

F(r+x)  =  F(x)  +  e  AxF(r) ,  V  r,  x  >  0, 

lii)  =  JIiL  ,  v  r,  x  >  o. 


and  thus 
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Hence  F(x)  =  c(l-eAx),  x  >  0,  and  f(x)  =  cAe  Ax  a.e.  on  (0,®).  In  view  of  the 
symmetry  of  the  distribution  of  X,  the  finite  constant  a  =  cA  may  be  taken 
positive.  Thus  (5.2)  is  shown. 

Conversely,  it  is  clear  that  the  process  (5.2)  satisfies  the  necessary  and 
sufficient  condition  for  being  left  weakly  Markov;  is  invertible,  in  fact  it  is  the 
stationary  solution  of  the  stochastic  differential  equation  X ’ ( t )  +  AX(t)  = 
aL1  (t) ,  where  L1  (t)  is  SaS  white  noise;  and  is  Markov,  as  follows  from  (s  <  t) 


(5.3) 


X(t)  =  e'A(t_sh(s)  +  e~A(t_u)dL(u) 


where  the  second  term  is  independent  of  {X(u),  u  <  s}. 

V'e  finally  show  that  the  process  (5.2)  is  the  SaS  Ornstein-Uhlenbeck  process 


v(t)  =  L(e“At>- 

(aA)1/a 


Indeed 


P(t)||“  =  aa/^e-°A(t-s)ds  = 

n-At  »a 

^ e  =  lx  ’ 

and  thus  ./{X(t)}  =  J^{Y(t)}.  From  (5.3),  using  the  independence  of  the  two  terms 
on  the  right  hand  side  we  have  for  s  <  t, 

8  (exp[irX(t)]  | X (s ) }  =  exp{ ire-A^t-s^X(s)  -  |ra|a  fl  eaA^t-u^du} 

=  exp{ire-A(t-s)X(s)  -  |r|aaQ'(l-^aA(t_s))/(aA)}. 

On  the  other  hand,  from  (2. 13)  with  H(t)  =  ae  At(aA)  ^a,  r( t)  =  eaAt,  we  have 

8  {e.xp[irY(t )]  |  Y(s) }  =  exp{  ire  A^  sW(s)-|r|aaa(  1-e  aA^  s)  )/(aA) ) . 


Therefore  J^{X(t)  | X (s) }  =  ./{Y(t)  | Y ( s ) }  for  all  s  <  t.  Since  both  X  and  Y  are 
Markov  it  follows  that  y(X)  =  y(Y).  It  is  clear  that  the  SaS  Ornstein-Uhlenbeck 
processes  Y  exhaust  the  class  of  all  SaS  Ornstein-Uhlenbeck  processes  (cf. 
(2.14)).  □ 

When  X  is  a  ScvS  moving  average  and  f  vanishes  on  the  positive  line,  we  say  that 
X  is  fully  anticipatory  and 

X(t)  =  /*  f(t-s)  dL(s) . 

A  fully  anticipatory  moving  average  is  called  invertible  if  sp{X(t),  t>r)  = 
sp{AL(t),  t>r}  for  all  r,  i.e.  the  increments  of  L  represent  the  backward 
innovations  of  X.  It  is  easily  seen  that  X(t)  is  fully  anticipatory  and  invertible 
(with  kernel  f(-))  iff  X  (-t)  is  nonanticipating  and  invertible  (with  kernel 
f (—  • ) )  •  Thus  the  only  fully  anticipatory,  invertible  SaS  moving  averages 
(1  <  a  <  2)  which  are  right  weakly  Markov  are  the  Markov  processes 

(5.4)  Y(t)  =  a'  /”  eV(t-s)dL(s),  <  t  <  «, 
where  0  <  a' .  A '  <  which  have  covariation  function 

Ry(0)  e_A'  t  >  0, 

(5.5)  Ry ( t )  =  «  I | 

Y  iRy(0)eAt,  t  <  0. 

Theorem  3.2  (iii)  implies  that  the  two  classes  of  Markov  processes  introduced  in 
this  section,  the  nonanticipating  and  the  fully  anticipatory  invertible  SoS 
moving  averages,  are  clearly  distinct  when  1  <  a  <  2. 
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6.  HARMONIZABLE  PROCESSES 

A  complex  harmonizable  S aS  process  has  a  harmonic  spectral  representation 

(6.1)  X(t)  -  /^eitudZ(u),  -oo  <  t  <  », 

where  Z  is  a  complex,  independently  scattered,  isotropic  SaS  measure  on 
with  p  a  finite  symmetric  measure.  For  every  complex-valued  function  g  e  Lft(//)  the 
r .  v.  / gdZ  is  complex  isotropic  SaS  with  if  exp{ SltzJgdZ}  =  exp{-| z|  a  /  |g|  ftd/z} ,  for 
complex  z.  Covariation  is  given  by 

Covf/gjdZ,  /g2dZ]  =  fg^^^dp 

for  all  gj,  g>2  €  L^(p),  where  z<ci>  =  j z | ^  ^z.  All  properties  of  covariation  and 
regression  used  here  in  the  real  case  are  also  valid  in  this  complex  isotropic  case 
(see  [4]) . 

The  harmonizable  process  (6. 1)  is  stationary  and  has  covariation  function 

(6.2)  R(t)  =  /^eitudp(u). 

Taking  its  real  part  would  provide  a  real  harmonizable  process,  but  it  is  more 
convenient  to  work  with  complex  quantities  when  dealing  with  Fourier  type 
representations.  Note  that  when  a=2  all  (continuous  in  probability)  stationary 
Gaussian  processes  are  harmonizable,  while  when  1  <  a  <  2,  the  harmonizable,  the 
nonanticipating  movirg  averages  and  the  sub-Gaussian  SaS  processes  form  disjoint 
classes  [6].  We  now  show  that  (nontrivial)  harmonizable  processes  cannot  be 
weakly  Markov. 

Theorem  6.1.  The  only  harmonizable  SaS  process  with  1  <  a  <  2  which  is  left  or  right 


weakly  Markov  is  the  constant  SaS  process. 


Proof .  Let  X  be  harmonizable  SaS  as  in  (6.1).  Assume  furthermore  that  X  is,  say, 
left  weakly  Markov.  Then  R(t)  =  R(0)e  for  all  t  >  0  and  some  0  <  A  <  ®.  The 
continuity  of  R  in  (6.2)  excludes  the  value  A  =  ®.  When  A  =  0,  X  is  a  constant  SaS 
process  with  /i  concentrated  at  0.  Assume  from  now  on  that  0  <  A  <  ®.  Since  R  is  an 
even  function  (cf .  (6.2))  it  follows  that 

R(t)  =  R(0)  e-A|t|  for  all  t, 


and  by  (6.2), 


du 


AR 


• 


tt(  A  +u  ) 

Since  X  is  right  weakly  Markov,  it  satisfies 

Cov[X(r)  -e"ArX(0),  Y]  =0,  Vr>0,  V  Y  6  sp{X(s) ,  s  <  0} 


Taking  Y  =  X(0)  +  X (— v )  =  ( 1+e  lvu)dZ(u)  with  v  >  0,  we  obtain  for  all  r,  v  >  0, 


0  =  F  (e  -e  )  ( l+e~  )  d/i(u) 
~00 

,»  (einl-e-Ar)(neivu)  du 


=  / 


[2( 1+cosvu)] 


1  a! 2  A2+u2 


Introduce  for  each  v  >  0  the  measure  (i  by 


d/iy(u) 

Hu 


(1+cosvu)^  Q//2  (A2+u2) 


2  /o 


Then  u  is  symmetric  and  finite  since  for  vu  5;  (2k+l ) ir,  1  +  cosvu  %  [vu— (2k+l ) tt]  /2 

V 

and  2(  1  — a/2 )  =  2  -  a  <  1 .  We  have  for  all  r,  v  >  0, 


/^xeinj(l+eivu)d/iv(u)  =  e_Ar/^(l+eivu)dMv(u). 


Let  f  (•)  be  the  characteristic  function  of  the  probability  measure  nv/nv( ER1) 
Since  nv  is  symmetric,  fy  is  real  and  even,  and  for  all  r,  v  >  0, 


fy(r)  +  fy(r+v)  =  e 


-A  t  1 


r  (l+cosvu)d/i  (u)  =  e  b(v). 


It  follows  that  for  all  k=l  ,2, . . . ,  and  v,  r  >  0, 


f v(r+2kv)  =  fy(r)  -eArb(v)  (l-e_Av+e~2Av- 

,  ,  -2kAv 

1+e 


_e-(2k-l)Avj 


and  thus 


limf  (r+2kv)  =  fy(r)  -e  Ar  b-^)-  ,  V  r,  v  >  0. 
k-oo  v  v  1+e  Av 


But  since  f  is  the  Fourier  transform  of  an  absolutely  continuous,  finite  measure 
by  the  Riemann-Lebesgue  lemma,  f(±oo)  =  0.  Hence  fy(r)  =  e  Arb(v)/(l+e  Av) 
V  r,  v  >  0,  and  by  the  symmetry  of  fy, 

fv(r)  =  e~A|r|,  V  r,  V  v  >0. 

v  1+e  Av 

The  uniqueness  of  Fourier  transform  now  implies  that  for  any  v  >  0, 


(1+cosvu; 


'(A4+u4)/iv(r)  1+e  tt(A%u4) 


— tj — 7j~,  for  almost  any  u 


and  thus 


(1+cosvu)1  =  — 7r(1+- - ,  ueK1, 


A  b (v)mv(Kx) 


which  is  a  contradiction,  as  the  left  hand  side  depends  on  u  while  the  right  hand 
side  does  not.  It  follows  that  X  with  0  <  A  <  <»  cannot  be  right  weakly  Markov  and 
similarly  it  cannot  be  left  weakly  Markov.  □ 

The  proof  of  Theorem  4.1  provides  an  example  of  a  stationary  SaS  process  with 


covariation  function  R(t)  =  R(0)e 


-A|t| 


(0  <  A  <  oo)  which  is  neither  right  nor  left 


weakly  Markov.  In  fact  this  is  a  harmonizable  process  with  representation 


(6.3) 


X(t)  =  [R(0)A/i]  /»  eitu  »  ■*,  dZ(u) 

(u%A  V' 


where  Z  is  complex  isotropic  SaS  motion. 
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7.  A  FAMILY  OF  ONE-SIDED  WEAKLY  MARKOV  PROCESSES 

Every  left  weakly  Markov  process  constructed  in  earlier  sections  is  also 
right  weakly  Markov,  and  vice  versa.  In  this  section  we  construct  non-Gaussian  SoS 
processes  which  are  left  weakly  Markov  processes  but  are  not  right  weakly  Markov, 
and  vice  versa.  Thus  for  non-Gaussian  SaS  processes  weak  Markovianness  is  not  a 
symmetric  property.  For  simplicity  we  consider  the  case  of  stationary  processes. 

Lemma  7.1.  Let  and  Xg  be  independent  left,  ( correspondingly  right)  weakly  Markov 
stationary  SaS  processes  with  covariation  functions  satisfying 

R;(t)  =  rj  e-3-!1! 

for  all  t  >  0  ( correspondingly  all  t  <  0) ,  i=l  ,2,  and  some  r^  a  >  0.  Then  for  any  real 
numbers  bj,  b2,  the  process 

X(t)  =  b^  Xj(t)  +  b2  X2(t) ,  —oo  <  t  <  oo, 

is  left  ( correspondingly  right)  weakly  Markov  stationary  SaS  process. 

Proof .  Only  weak  Markovianness  requires  proof.  We  assume  that  the  original 
processes  are  left  weakly  Markov.  For  any  r  >  0,  any  Uj,...,u^  <  0,  any  real 
C| , . . .  ,Cj.,  we  have,  using  the  independence  of  Xj  and  X2 , 

Cov[X(r)-e-arX(0),E^1ciX(Ui)] 

*  CovfbjfXjtr)  -  e-arX,(0))  *  b2{X2(r)  -  e“arX.2(0)}, 

bt^=l  CiXi(„i)  .b2^=iCiX2(ui)] 

=  |b,|0Cov[X1(r) -earX,(0),  sj[=1  C;  X,^)) 

+  |b2|aCov[X2(r)-e"arX2(0),  S^=1  X2(ui)]  =  0 


V'  *.  *.  V’  *.  v'  v" 


because  of  the  left  weak  Markovianness  of  and  X2.  Theorem  2.1  implies  that  X  is 
left  weakly  Markov.  □ 

Recall  that  we  have  introduced  three  families  of  weakly  Markov  stationary  SaS 
processes.  The  SaS  Ornstein-Uhlenbeck  processes  (2.13),  the  inverted  SaS 
Ornstein-Uhlenbeck  processes  (3.3),  and  the  sub-Ornstein-Uhlenbeck  processes  in 
Section  4. 

Fix  a  >  0.  Let  X1  be  SaS  Ornstein-Uhlenbeck  process  such  that 


Then  by  (2. 15)  we  have 


R1(t)  =  e  at,  t  >  0. 


tj(t)  =  e(Q^)at,  t  < 


Similarly,  let  X>2  be  inverted  SaS  Ornstein-Uhlenbeck  process  with 

R2(t)  =  e~at,  t>0. 

Then 

R2(t)  =  t  <  0. 

Finally,  let  Xg  be  sub-Ornstein-Uhlenbeck  process  with 


Then  (4.3)  implies  that 


R3(t)  =  e  at,  t  >  0. 


R3(t)  =  eat,  t  < 


The  three  processes  Xp  X2,  X^  are  assumed  to  be  independent.  Let  b^,  b2,  b^  be 
nonnegative  numbers,  at  most  one  of  which  is  equal  to  zero.  Define  a  new  stochastic 


process 


(7.1) 


X(t)  =  b^X^(t)  +  b2X2(t)  +  bgXg(t) ,  -00  <  t  < 


it 


By  Lemma  7.1,  X  is  left  weakly  Markov  stationary  SaS  process.  Using  the 
independence  of  X^,  X2 ,  X3  we  obtain 

R(t)  =  Cov[X(t) ,X(0)] 

=  bfCovp^t),  X1(0)]  +  bJCov[X2(t),X2(0)]  +  b"Cov[X3(t),X3(0)] 

=  b^R1(t)  +  b2  R2(t)  +  b3  R3(t) 


(b®  +  b2  +  b3)e  at  ,  t  >  0, 

+  bJe^/^+bJe111  ,  t<0, 


and  this  is  not  the  covariation  function  of  a  right  weakly  Markov  stationary  SaS 
process  (cf.  Section  2).  Therefore,  (7.1)  defines  a  whole  family  (with  parameters 
a,bpb2,b3)  of  left  weakly  Markov  stationary  SaS  processes  that  are  not  right 
weakly  Markov. 

It  is  clear  that  in  a  similar  way  we  can  define  a  family  of  right  weakly  Markov 
SaS  processes  which  are  not  left  weakly  Markov. 
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